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ABSTRACT 

Context. Molecular hydrogen is the main constituent of circumstellar disks and could be an important tracer for the evolution and 
structure of such disks. So far, H 2 has only been detected in a few disks and only through spectroscopic observations, resulting in a 
limited knowledge of the spatial distribution of the H 2 emitting gas. 

Aims. We report the detection of quiescent H 2 emission in a spatially resolved ring-like structure within 100 AU of T Tau N. We 
present evidence to show that the emission most likely arises from shocks in the atmosphere of a nearly face-on disk around T Tau N. 
Methods. Using high spatial resolution 3D spectroscopic K-band data, we trace the spatial distribution of several H2 NIR rovibrational 
lines in the vicinity of T Tau N. We examine the structure of the circumstellar material around the star through SED modeling. Then, 
we use models of shocks and UV+X-ray irradiation to reproduce the H2 line flux and line ratios in order to test how the H2 is excited. 
Results. We detect weak H 2 emission from the v=l-0 S(0), S(l), Q(l) lines and the v=2-l S(l) line in a ring-like structure around T 
Tau N between 071 (~ 15AU) and 077 (-100AU) from the star. The v=l-0 S(0) and v=2-l S(l) lines are detected only in the outer 
parts of the ring structure. Closer to the star, the strong continuum limits our sensitivity to these lines. The total flux of the v=l-0 S(l) 
line is 1.8 x 10~' 4 ergs s~'cirr 2 , similar to previous measurements of H 2 in circumstellar disks. The velocity of the H 2 emitting gas 
around T Tau N is consistent with the rest velocity of the star, and the H 2 does not seem to be part of a collimated outflow. Both shocks 
impinging on the surface of a disk and irradiation of a disk by UV-photons and X-rays from the central star are plausible candidates 
for the H 2 excitation mechanism. However, irradiation should not create a large degree of excitation at radii larger than 20 AU. Most 
likely the H 2 emission arises in the atmosphere of a flared disk with radius 85-100 AU and mass 0.005-0.5M o , where the gas is excited 
by shocks created when a wide-angle wind impinges on the disk. The H 2 emission could also originate from shock excitation in the 
cavity walls of an envelope, but this requires an unusually high velocity of the wide-angle wind from T Tau N. 

Key words. Stars: winds, outflows, circumstellar matter, Stars: emission-line, Stars: pre-main sequence, Infrared: stars 



1. Introduction 

The study of circumstellar disks around young stars is essential 
to understanding their evolution from gaseous disks to plane- 
tary systems. In this paper, we examine the spatial distribution 
of molecular hydrogen, the main constituent of disks around 
young stars. Disks have been observed in a wide range of wave- 
lengths ranging from optical to millimeter, although only a few 
studies have concentrated on the H2 component. Many inves- 
tigations have focused on the broad band spectral energy dis- 
tribution, which reflects the disk geometry and the structure of 
the dust content. Disks have also been observed more directly 
via optically thick dust lanes blocking the scattered light from 
young stars, and as near-infrared images of the scattered light o f 
the disk itself (e.g. lMcCabe et al.ll2002T:IWeinberger et al.ll2 002). 
Molecular line emission from species such as CO or HCO + is 
also used as a tracer for disks. The use of such tracers is, how- 
ever, subject to some uncertainty. Heavy element molecules may 
freeze out on dust grains, which likely settle to the midplane of 
the disk and/or get bound in larger rocks or planetesimals. Thus, 
molecules such as CO can become undetectable even if the disk 
still exists. 



1.1. Molecular Hydrogen in Disks 

Examining the H2 content in disks has many advantages. 
Hydrogen and helium are the last parts of the gas to be bound 
up when planets form, and will therefore remain in the disk after 
CO and dust have become undetectable. Observations of molec- 
ular hydrogen directly trace the gas mass of the disk without 
making assumptions about the dust-to-gas or CO-to-H2 ratios. 
Furthermore, molecular hydrogen will remain in the surface lay- 
ers of the disk when the dust settles to the midplane and is more 
directly accessible to incoming light than the dust and heavier 
elements. As a result, H2 may prove to be a better tracer for ex- 
ploring the evolution and structure of circumstellar disks, since 
it may be observable for a longer period of time. 

Direct observations of H2 in disks have been undertaken 
by several groups. We focus here on the observations of the 
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as well as fluoresce nt H2 in the UV I 
Herczeg et al.| [2006). Emission from the H2 v=l-0 S(l) line at 
2.1218/^m has been detected i n the disks of several T Tauri stars, 
classical as well as weak-line dBarv et al.ll2003l 120081: lltoh et al.l 
20031: IWeintraub et all 12001 iRamsav Howat & Greaves! 120071: 
Carmona et al. 120071) . These detections are made through 
longslit spectroscopic observations, and they do not reveal much 
about the spatial distribution of the molecular hydrogen, beyond 
indicating that the emitting gas is located within 50 AU of the 
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central star. IChen et all d!998l) presented images of H2 v=l-0 
S(l) emission from photoevaporating disks in Orion and showed 
that the emission arises on the disk surface. In this case, the disks 
were externally irradiated, and the H2 emission was found from 
a region ~ 200AU in size. 

In this paper, we present spatially resolved images of H2 
emission from a ring around T Tau N obtained with the integral 
field spectrograph SINFONI on the ESO-VLT. The presence of 
H2 emission in the T Tau system has been known for decades, 
but this is the first time that the weak emission within 100 AU of 
T Tau N has been resolved and analyzed. 

1.2. T Tau 

T Tau is a triple s tar system with an age of ~ lMyr 
dWhite & Ghezl [200 lb . The binary component T Tau S, con- 
sisting of T Tau Sa and T Tau Sb (separation ~071), is cur- 
rently ~077 south of T Tau N. Al l three stars are activ ely ac- 
creting and believed to host disks dDuchene et al.ll2005b . T Tau 
S shows heavy extinctio n (Ay = 15), which is attributed to 
a circumbinary structure (Duchene et al. 2005). Another possi- 
bility is that T Tau S is obscured by the disk around T Tau N 
(iHoeerheiide et al.lll997t iBeck et al.ll200lh 

T Tau N is a ~2M star dWhite & Ghezl200ll) and is believed 
to have a disk that is seen nearly face-on ( Akeson et al.| [l998). 
Based on photome tric periodicity and assume d stellar radius, 
iHerbstet al] d 19971) derive an inclination of 19°. IStapelfeldt et all 
( 19981) suggest an outflow and disk with the axis at position an- 
gle 300° and with inclination of ~ 4 5° in order to expla in the 
morphology of scattered optical light. Akeson et al. (2002) find 
the inclination to be 20 - 40° from SED fitting. 

This paper is organized as follows. In Sect. [2] we describe 
the observations and data reduction. In Sect. [3] we present the 
spatial distribution of molecular hydrogen around T Tau N and 
the velocity distribution of the gas. Section [4] discusses the ge- 
ometry of the star-disk-envelope system and Sect. examines 
the H2 excitation mechanism. In Sect. [6] we consider the possi- 
ble implications for T Tau S, and finally, we draw conclusions in 
Sect. E| 



2. Observations 

T Tau was observed with the ESO-VLT as part of the SINFONI 
science verification program on the nights of 2004 October 30th 
and November 2nd . SINFONI is a near-infrared (NIR) integral 
field spectrograph work ing in combination with adaptive optics 
(Eisenhauer et al. 2003). Observations of the region around the 
T Tau triple star system were obtained in the K-band using the 
3'.'2 field of view optics (100 mas pixel scale) centered on the 
northern component. T Tau N (my=9.6) itself was used as the 
guide star, producing diffraction limited spatial resolution. The 
2D image on the sky was sliced into 32 slitlets which were then 
dispersed onto a 2k x 2k detector. The spectrograph provides a 
spectral resolution of 4000 in the K-band. The observations were 
carried out using a five-point nodding pattern with individual ex- 
posure times of 3 seconds and 20 co-adds and a total integration 
time on source of 30 minutes. The nodding pattern was that of a 
box with one arcsecond width centered on T Tau N. The result- 
ing mosaic has a field of view of approximately 4 ". Sky frames 
with the same exposure times were obtained within the nodding 
cycle. 

Data reduction and reconstruction of the 3D cubes were car- 
ried out using the SINFONI pipeline (version 1 .3.0) provided by 



ESO. The 2D raw frames were corrected for sky background, flat 
field effects and optical distortions. Bad pixels and cosmic rays 
were identified and the frames were calibrated in wavelength. 
Then, the 3D cubes were constructed using calibration data of 
the positions and distances of the slitlets on the detector. The 
cubes within the nodding cycle were aligned spatially and coad- 
ded plane by plane to create the final mosaic. Since the total ex- 
posure time is less at the outer regions of the mosaic than in the 
centre we scaled the flux at all spatial points to an exposure time 
of 3 seconds. The final 3D cube stores the spatial information 
in the x- and y-directions and the spectral information along the 
z-direction. To improve the signal-to-noise ratio, each spectral 
plane was smoothed with a 3 by 3 boxcar in the spatial domain. 

The B9 standard star Hip025657 was observed under the 
same conditions and similar airmass as T Tau and with the same 
instrumental setup, in order to correct for atmospheric absorp- 
tion. The spectrum was extracted after the data had been re- 
duced following the same recipe as for T Tau. The spectrum 
of Hip025657 is featureless except for Bry in absorption. We 
removed this feature and replaced it by a linear fit to the sur- 
rounding continuum. Subsequently, the spectrum was divided by 
a blackbody function of T= 11000 K and normalized. Dividing 
each spectrum of the science cube by the corrected standard star 
spectrum removed telluric absorption features in the T Tau spec- 
tra very effectively. 

Flux calibration was also performed using Hip025657 
(m/(-=7.443). The conversion factor between counts s _1 and 



ergs s 1 cm 2 fim 1 sr 1 was found by dividing the K-band 
flux of the star (4.1 xl0~ 7 ergs s crrr 2 jurrT 1 



xl0 -»u/2.5 

dCampins et alj [19851) ) by the mean counts per second of the 
standard star spectrum within 1.94 - 2.45/^m and dividing by 
the pixel area in steradians. 



3. Results 

We show the spatial distribution of emission in the H2 v=l-0 
S(l) rovibrational line at 2.12jt/m in Fig.Ql The image is domi- 
nated by strong H2 emission south of T Tau N, close to T Tau S . 
Thi s emission feature was also detected by Herbst et al. (2007) 
and lBeck et all d2008l) . The origin of the strong emission south 
of T Tau N is belie ved to be outflows from one or more of the 
stellar components (Herbst et al. 2007) and will be the subject 
of a subsequent paper (Gustafsson et al. in preparation). In this 
paper, we focus on the weaker H2 emission found very close to 
T Tau N. The weak feature is seen to extend all around the star 
in a ring-like structure. The morphology of the H2 emission in 
our data is fully c onsistent with the map recently published by 
IBeck et al] d2008l) . They also detected the weak ring-like struc- 
ture although they did not mention it. Their data were obtained 
October 2005, which indicates that the emission feature is stable 
on at least a 1-year timescale. 

A close-up of the immediate surroundings of T Tau N ap- 
pears in Fig. [2] H2 emission is found as close as 071 arcsec (~ 
15 AU, assuming a distance of 140pc) and is seen to extend out 
to ~ 077 (100 AU) from the star. We do not detect H2 emission 
above the noise closer to the star than 07 1 . Molecular hydrogen 
does not appear in the spectrum of T Tau N (Fig. [3} and is an 
exclusively extended phenomenon. 

Emission from other H2 lines than the 1-0 S(l) transition has 
also been detected (Fig. [3]). The Q(l) line at 2.406//m is found 
to have roughly the same spatial distribution as the S(l) line. 
Other Q-branch lines are present as well, but the correction for 
atmospheric absorption is challenging in this spectral and spatial 
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Table 1. Total flux of H2 lines within mask. : 3<x upper limit. 
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Fig. 1. H2 v=l-0 S(l) emission in the T Tauri system. The color 
scheme indicate the flux level in ergs s _I cirT 2 sr -1 . The posi- 
tions of T Tau N and T Tau S are marked with white circles. The 
white box outlines the close-up view of the ring-like structure 
around T Tau N shown in Fig. [2] 
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region and may introduce errors. The S(0) (2.2223//m) also ap- 
pears in the spectrum. This line is, however, weak and the power- 
ful continuum emission from T Tau N limits the line detection at 
a level of 3<x above the noise level to pixels located in the outer 
ring structure at r> 70 AU. At larger radii, the S(0) emission 
shows a similar spatial distribution to the S(l) line. The v=2-l 
S(l) line (2.2447/im) is tentatively detected at a 2<x level at a few 
locations in the outer region of the ring structure at r> 80 AU. 

The total flux of the H2 v=l-0 S(l) in the structure is found 
by summing all light within a mask, the extent of which is shown 
in Fig. [2] We have chosen a rather conservative mask in order to 
avoid confusion with emission features that may have a different 
origin. Thus, the mask excludes regions where the emission is 
weaker than 4 x 10~ 4 ergs s" 1 cm" 2 sr~' as well as the strong 
emission region south-west of T Tau N which is most likely 
cause d by an outflow from one of the T Tau S stars dHerbst et alj 
120071) . The total flux within the mask is 1.8 ± 0.3 x 10~ 14 ergs 
s _1 cm" 2 . The uncertainty is the formal uncertainty on the total 
flux calculated using the flux uncertainty in the pixels included 
in the sum. This uncertainty does not take into account that the 
estimated flux depends on the chosen mask and that the exact 
shape and extent of the ring-like structure is difficult to quantify 
because of other H2 features nearby. We estimate that the un- 
certainty due to the mask may amount to 30%. The flux of the 
v=l-0 Q(l) line and the upper limits to the flux of other H2 lines 
appear in Table [TJ 

The total flux in the v=l-0 S(l) line is similar to the amount 
of H2 emission at the same radial velocity as the star detected in 
the circumstellar environment of other T Tauri stars, where it is 
believed to originate from disks within 100 AU (Bar vet al.l200l 



Fig. 3. K-band spectrum of T Tau N and of a pixel within the 
mask in Fig. [2] located (0775,0) with respect to T Tau N. Bry 
emission and photospheric absorption lines are seen in T Tau N. 
In the extended ring-like structure, only H2 lines are detected. 
The Bry line in the lower spectrum is contamination from the 
PSF of T Tau N and is not due to extended emission of Bry. 

I2008h I Weintraubetai] 120051: iRamsav Howat & Greaves! 120071) . 
The H2 line flux previously measured in disks ranges from 7 x 
10~ 16 - 1.5 x 10~ 14 ergs s" 1 cm" 2 . 

In our data, we have full access to the spectral distribution 
of the emitting gas. Although SINFONI only has a spectral res- 
olution of ~ 75 kms -1 in the K-band, it is possible to deter- 
mine the peak position of the lines with much higher accuracy 
through line fitting. We have derived the radial velocity corre- 
sponding to every H2 emitting position (see Fig |2]i by fitting a 
Gaussian profile to the unresolved line profiles on a pixel by 
pixel basis. The velocities in Fig.|2]have been corrected for the 
Earth's motion toward T Tau at the time of observation and are 
quoted with respec t to the heliocentric v elocity of T Tau N of 
19.1 ± 1.2kms _1 dHartmann et al.l IT986I) . The velocity map in 
Fig. [2] shows small velocity variations between -lOkms -1 and 
+ 10kms _1 with respect to the intrinsic velocity of T Tau N 
within the ring structure. Th ese velocities are consistent with the 
data from Bec k et all d2008l) who also measure velocities close to 
the systemic velocity at this location (see their figure 1 1). There 
is a tendency for a radial gradient with the velocities being posi- 
tive (5-lOkms -1 ) close to the star and negative (-5kms _1 ) further 
out. However, the uncertainty in the derived velocities are larger 
close to the star than further out due to the increased continuum 
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Fig. 2. Top: H2 v=l-0 S(l) emission surrounding T Tau N. The color scheme indicate the flux level in ergs s _I cirT 2 sr~' . Positions where no IT 
emission is detected are shown in black and the position of T Tau N is marked by a white circle. The mask outlining the region used for further 
analysis is indicated by the black line. The strong emission region south-west of T Tau N is excluded because it is most likely related to outflows 
from T Tau S (see text). Centre: The derived heliocentric radial velocity of the H2 emission relative to the rest velocity of T Tau N. Bottom: 
Uncertainties in the radial velocities. 
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Fig. 4. H2 v=l-0 S(l) line profile integrated over the ring struc- 
ture (diamonds with lcr errorbars). Fitted Lorentzian (dashed 
line). 



emission in the inner region (Fig. [2j. This makes any conclusion 
on the radial variation of velocities rather uncertain. There is no 
evidence of Keplerian rotation of the disk. However, if the incli- 
nation of the disk is the same as the star itself, ~ 20°, the radial 
velocity component of Keplerian rotation around a 2M star is 
only ~ 5kms _1 at 10AU and ~ 2kms~' at 100AU. Such small 
velocity differences within the disk would be difficult to detect 
with the present data. 

In order to improve the signal-to-noise ratio, we constructed 
a global H2 profile of the ring-like structure by adding all spec- 
tral profiles of H2 v=l-0 S(l) emitting positions within the mask 
in Fig. [2] This also allows a direct comparison with previous spa- 
tially unresolved measurements of H2 in the circumstellar envi- 
ronment of T Tauri stars. The global line profile appears in Fig. [4] 
together with a Lorentzian fit. The Lorentzian fitting function 
provides the best match to the instrumental profile of SINFONI 
which dominates the unresolved H2 profile. The profile is seen to 
peak close to the rest velocity of T Tau N. From the Lorentzian 
fit we find the peak velocity to be -2.5 + 2.1kms~ 1 (lcr un- 
certainty). Considering the uncertainty in the rest velocity of 
T Tau N of 1.2kms^ dHartmann et alj|1986l) . the velocity of 
the H2 emission is consistent with the rest velocity of the star 
within the errors. The same was foun d to be true of the H2 
emission from disks around other stars dBarv et al.ll2003l [2008 ; 
iRamsav Howat & G reaves 20071 : ICarmona et al.ll2007l) . 

4. Disk or outflow? 

We now examine the origin of the H2 emission around T Tau 
N. The ring-like shape of the emission can be created by several 
scenarios. We consider the following possibilities: 1) an outflow, 
2) shocks created by a wide-angled wind hitting the disk or cav- 
ity walls , 3) a photo-evaporating wind from a disk viewed almost 
face-on dHollenbach et al.lll994tlFont et al.ll2004l) , 4) a photodis- 
socia tion region, or 5) UV/X-ray heating of a disk ( Nomura et al. 
120071) . In this section, we consider the geometry of the system 
and in Sect. [5] we examine what excitation mechanisms can re- 
produce the H2 line flux and ratios. 

Outflows from T Tauri stars are known to be complex, 
often showing two velocity components with different spatial 
characteristics. Observations of both forbidden emission lines 



a 



b 



c 




Fig. 5. Excitation scenarios in the T Tau N-disk-envelope sys- 
tem, a): a wide-angle outflow hits the walls of a cleared out 
cavity in the surrounding envelope, b): a wide-angle outflow 
impinges on a flared disk and creates shocks, c): UV + X-ray 
irradiation from T Tau N heats a flared disk. 



dBacciotti et al] |2000|) and molecular hydrogen ( Takam Tet al] 
2006) have shown that outflows typically consist of a collimated 
high velocity (60-200kms _I ) jet and a less collimated low ve- 
locity (0-30kms _1 ) component, i.e. a wide-angled wind. At first 
glance, the small line-of-sight velocities measured in the ring 
argue against a collimated jet. Although a jet viewed almost 
pole-on could create a circular emission feature, the velocities 
associated with such a jet are much higher than observed here. 
It is, however, possible that we see shocks from a wide-angle 
wind interacting with either the outer walls of a bi-conical cavity 
cleared out by the outflow or a flared disk (Fig. [5^ or[5}3 ). The 
existence of a cavity in an enve lope around T Tau N was sug- 
ge sted by Momose et al] d 19961) and the reflected light images 
of ISta pelfeldt et al] d!998l) . In this picture, a wide-angle wind 
creates oblique shocks when interacting with the molecular en- 
vironment. Oblique shocks result in low shock velocities, since 
it is only the normal component of the velocity with respect to 
the gas that is thermalized. Furthermore, the measured veloci- 
ties will be even lower, since only the line-of-sight velocity is 
detected here. 
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Table 2. SED model parameters 



envelope accretion rate 


10"'Mo/yr 


cavity opening angle 


47° 


Disk mass 


0.15 M s 


Disk outer radius 


85 AU 


Disk inner radius 


0.3 AU 


flaring parameter, [} 


1.14 


Disk accretion rate 


10- 6 M o /yr 


Disk inclination 


< 65° 



Another possibility is that the emission is linked to irradia- 
tion of a nearly face-on disk around T Tau N (Fig[5J;). We pro- 
ceed by investigating the circumstellar environment of T Tau N 
through spectral energy distribution (SED) modeling to find ev- 
idence of the size of the disk and envelo pe. We use the precom - 
puted grid of radiative transfer models of Robitai lle et al.l (|2007). 
Their models include contributions from a circumstellar accre- 
tion disk, an infalling envelope and an outflow cavity. The mod- 
els span a wide range of parameter space. The SED is calculated 
for each set of model parameters for comparison with observa- 
tional data. Using fluxes of T Tau N from the literature, we fit 
the optical and infrared data with the models in order to find the 
range of physical parameters that fit the SED best. 

The SED appears i n Fig. [6] Th e observations include 

Kenyon & Hartmann 



(fl995b: 



(119971 T2 007): 
Hog erheiide et all dl997b; 



those of [Herbst et al 



Weaver & Jones (1992); 



lAndrews & WiUiamsl(l2005h : lBeckwith & SargenJ(ll99lh . T Tau 
N shows brightne ss varitations of 0.2 mag in the K and L bands 
(Bec k et ai1l2004b . which can explain the spread in flux values 
that the various authors find. We have included measurements 
from different epochs at optical and infrared wavelengths, in or- 
der to average out the effect of this variability. In the model fit, 
we have only used data points at A < 21/im, for which the T Tau 
N component is resolved. At longer wavelengths, the north and 
south components o f the T Tau triple system are unresolved and 
T Tau S dominates dGhez et alj|199lh . We assume that Ay = 1.5 
dWhite & Ghezll200ll) . 

The Robitaille grid resulted in a wide range of parame- 
ters that produced a reasonable fit to the SED. In order to fur- 
ther constrain the parameter space, we used the spectral infor- 
mation on T Tau N, a ~ lMyr o ld K0 star with M~ 2M , 
T eff = 5200K dWhite & Ghezll2001l) . Thus, we chose to consider 
only those models in which the stellar parameters fall within 
the following ranges, 1.5 <M/M <2.5, 4000K<T eff <6000K, 
0.1Myr<age<2Myr. The best fit model which satisfies these cri- 
teria appears in Fig. [6] For the 50 models that best fit the data, 
only four models satisfy the stellar restrictions of T Tau N. These 
four models have the same values for all parameters except the 
inclination angle of the disk. The model parameters are listed in 
Table [2] The inclination is poorly determined, but is found to be 
less than 65°. 

SEDs only contain a limited amount of information about the 
distribution of circumstellar material, and it is easy to overinter- 
pret the results of model fitting. Nevertheless, our results show 
that the T Tau N system is consistent with an accreting star sur- 
rounded by a disk and torus-like envelope with an opening angle 
of ~ 45°. The disk is likely to have an outer radius of 85 AU con- 
sistent with the detection of H2 at a radius of ~ 100AU. It also 
shows a large degree of flaring. The inclination of the disk cannot 
be well constrained by the SED-model ing but is consisten t with 
a low inclination of ~ 19°, as found by Her bst et"aT1 dl997[). The 
model is also consistent with the suggestion oflStapelfeld t'eTal] 
(1998) that an outflow at PA=300° and inclination of 45° has 



blown out a cavity. Note that the disk accretion rate is higher 
than the observed accretion rat e of (1.4 - 5.9) X 10~ 8 M G yr _1 
based on the Bry line strength dBeck et al.ll2.004t >. This is most 
likely due to inconsistencies in the radiative transfer models used 
in the Robitaille models, which systematically overestimate the 
accret ion rate. This effect was already noted in Robitaill e et al.l 
(2007). Thus, the accretion rate inferred from the SED modeling 
should be considered with some caution. 

Assuming o ptically thin emis sion, the mass of H2 v=l-0 S(l) 
emitting gas is (Bary et al. 20031) 



M(H 2 ) 



f=1-05(1) 



i) = 1.76 X io-x> 4 f^ DZ = 5 x 10- 10 M e , (1) 



where Fu ne is the v=l-0 S(l) flux, D the distance in pc, E v j the 
upper level energ y and Au ne the Eins tein A coefficient of the tran- 
sition. Following iBary et all d2003l) . we estimate the total disk 
mass by multiplying the v=l-0 S(l) H2 mass by a scaling factor 
of 10 7 - 10 9 . This scaling factor was found by comparing the H2 
mass to masses derived from CO and (sub)millimeter continuum 
in disks around four T Tauri stars (GG Tau, LkCa 15, DoAr 21, 
TW Hya). This gives a disk mass of 0.005 - 0.5M Q . Note that 
the most likely disk mass of 0.1 5M Q estimated from the SED 
modeling falls within this range. Analysis of HCO + and mil- 
limeter continuum data in dicate that the disk mass around T Tau 
Nis 1 0' 3 - 4x 10~ 2 M Q dHogerheiide et al.|[l997l;lAkeson et al.l 
1998), suggesting that the actual disk mass is at the lower end of 
the simple mass estimate based on the H2 mass. 

These arguments suggest that T Tau N has a flared disk with 
radius 85-100 AU and mass 0.005 - O.5M . We propose that 
the observed H2 emission comes from the gaseous disk atmo- 
sphere or alternatively from the torus-like envelope that is also 
present around T Tau N according to the SED. The reason that 
T Tau N is not at the centre of the H2 emission (Fig. |2]i could 
be due to the fact that we do not see the disk-envelope system 
exactly face-on. An inclination of 20-40° would cause emission 
from a circular structure to appear elliptical. Assuming that we 
only see emission from the surface facing toward us projection 
effects would furthermore make the nearest side of a disk appear 
narrower than the farther side. This would shift the projected po- 
sition of the central star away from the centre of emission. A 
detailed modeling of the spatial distribution of the emission is, 
however, outside the scope of this paper. 

If the value of 85-100 AU reflects the true size of the disk, 
why has the disk not yet been detected in scattered light? The 
simplest explanation is that the disk is small (r < 077) and the 
flux contrast between the star and disk is large, making detec- 
tion difficult. Infrared observations of scattered light from nearly 
face-on disks (circumstellar and circumbinary) around other T 
Tauri stars have shown th at the total flux from a disk is only 1-2 
perce n t of the stellar flux dMcCabe et al . 2002; Weinb erger et al.l 
12002b . IStapelfeldt etail (11998b and iMavama et al.l d2006h ob- 
served reflected light in T Tau at optical and infrared wave- 
lengths, respectivel y, but did not find evid ence of the disk around 
T Tau N. However, Ma vama et al.l (12006b used a coronographic 
mask with a diameter of ~075 -076, which is only marginally 
smaller than the size of the disk found in this paper. The disk 
would therefore in any case be difficult to see in their da t a. The 
detection of the disk from the data of IStapelfeldt et al.l ([1998) 
could be compromised by artefacts arising from the PSF sub- 
traction. 
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Fig. 6. Spec tr al energy distribution o f T Ta u N and model cur v es. Th e data points inclu d e thos e of this work and Herbst et"ai] 
d!997l 120071): iKenvon & Hartmannl (Il995l) : iHogerheiide etail dl997h : IWeaver & Jonesl (Il992l) : lAndrews & Williamsl d200ll) : 
Beckwith & Sargent (1991). Plus-signs indicate resolved data of T Tau N, diamonds indicate where T Tau S and T Tau N are 
unresolved and the flux of T Tau S dominates. The K-band flux from this work is marked with a triangle. The model is found from 
fitting to data at A < 21/im where T Tau N and T Tau S are resolved. Solid line: total flux, dotted line: stellar flux, dashed line: disk 
flux. 



5. Excitation mechanism 

After studying the location of the circumstellar material and 
the H2 emitting gas, we now turn to the H2 excitation mecha- 
nism. The near-infrared rovibrational lines of H2 can be excited 
through different processes. The two main mechanisms are (i) 
shocks (Fig(5^ and [5J3) and (ii) UV + X-ray radiation (Fig|5j;). 
In the former case, the H2 molecules are thermally excited by 
the passing shock wave and the H2 spectrum is characterized by 
a single excitation temperature. In the latter case, the UV and X- 
ray radiation may contribute both to the heating of the gas and to 
electronic excitation of H2. The near-infrared H2 emission may 
therefore contain a contribution from both thermal excitation 
as well as non-thermal radiative decay from excited electronic 
states (fluorescence). These two excitation mechanism can be 
distinguished through both the line strengths and the line ratios 
of the observed H2. Traditionally, the ratios /,.=i-05(i)/A=i-os(0) 
and /r=i-os(i)/A'=2-i5(i) are used. These ratios also provide a di- 
agnostic for whether the thermal or non-thermal contribution 
dominates, since higher vibrational bands, v > 2, are more 
densely populated by fluorescence than by thermal excitation. 
In the following, we will concentrate on the two excitation pro- 
cesses (i) shocks and (ii) UV + X-ray irradiation 



5.1. Shocks 

In the first process, shocks, a super- Alfvenic shock wave rapidly 
heats the gas to temperatures of ~ 1000 - 3000K. The temper- 
ature reached depends on the velocity of the shock, as well as 
the pre-shock conditions in the ambient medium, such as den- 
sity, magnetic field, chemistry etc. The /,.=i-0S(i)/A>=2-is(i) ratio 
depends on the temperature, since the populations are assumed 



to be in local thermodynamic equilibrium (LTE). The higher the 
temperature, the larger the relative population in higher vibra- 
tional bands. We estimate the range of the I v =i-os(i)/Iv=2-iS(i) 
ratio that shocks are able to produce considering the temperature 
range 1000 -3000K reached during the passing of the shock. For 
an optically thin transition, the line intensity is 



Au ne hc 

lline - — 7-jy v J, 

AnA 



(2) 



where N v j is the population in the upper level denoted by the 
v and J quantum numbers, A the wavelength, and Au„ e is the 
Einstein coefficient of the transition. We ignore extinction ef- 
fects, since the differential extinction between the H2 lines is 
small. In LTE, the level population is given by 



N v _ 



N„ 



V..I 



Z(T) 



g vJ exp(-E vJ /T), 



(3) 



where N lo , is the total column density, Z(T) is the partition func- 
tion, g v j is the degeneracy of the level, E v j is the level energy 
expressed in K, and T is the excitation temperature. The degen- 
eracy depends only on the rotational state, J. Thus, g\$ = #2,3- 
We find: 



^•=i-05(i) A,, = i_os(i) 2.24/im (E23-. 

= exp\ 

/v=2-is(i) A,. =2 -i5(i) 2.12yuw \ T 

= 0.736^1^^1, 



•1,3 



(4) 



which gives a value of 198 at 1000 K and 5 at 3000K. Similarly, 
we find / v =i-05(i)Mv=i-os(0) - 3 - 5 for T~ 1000 - 3000K. 

From these estimates, it is clear that a range of shock tem- 
peratures can reproduce the line ratios observed around T Tau 
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N (TableQ]). However, since the observations give us only lower 
limits to the line ratios, we are not able to constrain the shock 
temperature and thus the underlying physical parameters very 
well. We therefore turn to the total line flux. Shock models show 
that the H2 v=l-0 S(l) surface brightness in Fig. [2] can be pro- 
duced in different density regions of the circumstellar material. 
The required density is tightly correlated with the impact veloc- 
ity of the shock front. This means that the H2 emission can arise 
in environments with pre-shock density of n# ~ 1 x 10 4 crrT 3 if 
the shock velocity is ~ 40 - 50kms _1 while a shock front prop- 
agating at only ~ lOkms -1 is sufficient to excite the H2 if the 
pre-shock density is n# ~ 1 x 10 7 cirT 3 (L. E. Rristensen, private 
communication). 

From models of T Tauri stars with disk and envelope masses 
resembling those of T Tau N, it is evident that the high density 
regime of iih ~ 1 x 10 7 crrr 3 is associated with the upper layers of 
the disk while densities of n# ~ 1 x 10 4 cirT 3 are found in the en- 
velope and in a thin transition zone between the disk atmosphere 
and the outflow cav ity, i.e. the cavity walls dWhitnev et al.ll2003"t 
Crap si et ai1l20 08). Independent of where the excitation takes 
place, the outflow or wind must impact the gas at a small angle, 
unless the disk is highly flared. Thus, the flow velocity must be 
much higher than the shock velocity, since it is only the velocity 
component perpendicular to the shock surface that contributes to 
the shock. Assuming an impact angle of 10°-20°, a flow velocity 
of 120-230kms _I is required to create a shock of 40kms _1 . Such 
velocities are found in the collimated jets from T Tauri stars. 
However, the SED indicates that the cavity opening angle is large 
(~ 47°, Table |2]i and thus it is difficult to see how a collimated 
jet should be able to hit the cavity walls. On the other hand, a 
shock velocity of lOkms -1 can be attained with a flow velocity 
of 30-60kms _1 and an impact angle of 10°-20°. Such velocities 
may be found in a wide-angled low velocity wind charac teristic 
of T Tauri stars dBacciotti et al. 2000; Taka mi et alj |2006). 

In summary, if shocks are the primary cause of the H2 exci- 
tation, it seems more likely that they arise from the interaction 
of a wide-angled wind with the upper layers of a flared disk than 
with cavity walls carved in a circumstellar envelope. We cannot 
exclude that the H2 excitation arises from the cavity walls, but 
this scenario seems to require that the high velocity jet in T Tau 
N be less collimated than in other sources. 

5.2. Irradiation 

We now turn to excitation process (ii), irradiation by UV-photons 
and X-rays (FiglSj;). In the i r mod els of fluorescent excitation of 
H jBlack & va n DishoecS dl987l) find/ v= i_os(i)//v=i-os(0) < 2.7 
and /v=i-05(i)//v=2-is(i) < 2.0. This is incompatible with our 
data (Table[TJ and shows that fluorescence alone cannot explain 
the observations. A thermal component is nece ssary in order 
to reproduce the line ratios. iLe Petit et al.l d2006l) have shown 
that the ratio /v=i-as(i)/^v=2-is(i) can be much larger than 2 in 
a photon-dominated region (PDR) if n# > 10 5 cm 3 and the 
incident far-ultraviolet radiation field is stronger than ~ 10 4 
times the average in the interstellar medium. In such high den- 
sity regions, collisions of H2 in vibrationally excited states re- 
sulting from fluorescence with atomic H tend to thermalize the 
rovibrational states. These conditions may very well apply to a 
dense circumstellar disk irradiated by the central star. The bright- 
ness of the v=l-0 S(l) line in the model is found to be larger 
than ~ 5 x 10~ 4 erg s s _1 cirT 2 sr _1 in a PDR viewed face on 
dLe Petit et al.ll200a) . which is consistent with our data. The H2 
brightness around T Tau N is larger than this value out to dis- 



tances of ~ 80AU (Fig. [2]). The H2 emission may thus be caused 
by UV-irradiation from T Tau N creating a PDR at the dense sur- 
face of the disk. Such a PDR would create fluorescent H2 lines 
in the UV and infrared H2 lines composed o f contributio n s from 
both fluorescence and thermal excitation. IWalter et al.l (2003) 
found extended fluorescent H2 emission i n the UV around T Ta u 
N. This seems to support our conclusion. iHerczeg et al.l d2006). 
however, did not find fluorescent H2 emission in the UV extend- 
ing more than O'.'l from T Tau N.[Saucedo et al. (2003) reported 
fluorescent H2 emission North-East and South- West of T Tau N 
at a distance of < 50 AU, which is consistent with pumping by 
stellar Lya emission. They do, however, find that the fluorescent 
line needs to be heated before being pumped and suggest that an 
outflow must be the heating mechanism. 

iNomura et al.l d2007l) constructed models of H2 emission 
from a disk irradiated by both X-rays and UV photons from a 
central T Tauri star. For a model with a 0.5M o , T e ff - 4000K 
central star with X-ray luminosity of Lx ~ 10 3() ergs s _1 and 
UV excess like TW Hydrae, they find v=l-0 S(l) line fluxes 
of 0.1 - 20 x 10~ 15 ergs s _1 cirT 2 , depending on the size 
of the dust grains. When the grain size increases the high- 
temperature region in the disk shrinks and the line flux de- 
creases. The /,.=i-05(i)/A=2-i5(i) ratio is always larger than 4 
and can be as high as 50 if UV radiation dominates. In units 
of 10~ 14 ergs s _1 cirT 2 , the estimated flux in the v=l-0 S(l) 
line for a model using small dust grains is found to be ~ 0.1 
when only X-ray irradiation is considered, -1.3 with only UV- 
irradiatio n and ~ 2.1 when b oth X-rays and UV-irradiation are 
included dNomura e t al. 2007). The estimated total flux from the 
disk around T Tau N is 1.8 on this scale, suggesting that X- 
rays alone are not sufficient to produce the observed H2 emis- 
sion. UV-photons seem to dominate the excitation but a combi- 
nation of X-rays and UV-photons may be necessary in order to 
explain the H2 flux as emission from a disk around T Tau N. 
Note that T Tau has an X-ray luminosity of Lx = 3.5 x 10 31 ergs 
s" 1 dGiidel et alj2007l) and UV-luminositv of Lnv = 3xl0 33 ergs 
s" 1 , c orresponding to a large UV-excess of 0.7L o (ICalvet et al.l 
12004 . 

It seems likely, then, that the UV-irradiation from T Tau N 
itself, with a possible X-ray contribution, is strong enough to 
produce the to tal H2 flux we see around T Tau N. However, in 
their models. INomura et al. I (12007b find that most H2 emission in 
the v=l-0 S(l) line is emitted at a radius of 20AU and that the 
emissivity decreases with increasing distance from the star. The 
presence of H2 emission extending to ~ 100AU around T Tau 
may be difficult to explain with models of irradiation. A possible 
way to solve this is if the disk is strongly flared. A large degree 
of flaring will allow more UV-photons to reach the outer parts of 
the disk, thus increasing the excitation in the outer regions. 

5.2.1. Photoevaporation 

If the disk is irradiated by UV photons and X-rays, 
a disk-win d may be powered by photoevaporat i on of 
the disk fljo llenbach etaf] Il994t Uohnstone et all I1998I; 
IStorzer & Hollenbachlll999l: iFont et ai1l2004l) . The high-energy 
radiation (EUV + X-rays) ionizes hydrogen at the disk surface 
and heats the gas to 10 4 K. EUV photons cannot penetrate 
the ionization fron t and X-rays heat o nly the inner disk and 
the surface layer (No mura et alj 120071) . FUV photons, which 
are not absorbed by atomic hydrogen but mainly by dust, 
penetrate much deeper and reach the disk surface. Here, they 
dissociate molecular hydrogen and heat the neutral gas to 
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about 400-4000 K. If the thermal velocity of the neutral gas 
exceeds the escape velocity of the disk surface, the gas flows 
outward. The flow is initially cylindric al but is reorien ted into a 
spherical flow by pressure gradients dFont et al.ll2004l) . Deeper 
in this photon-dominated region, there is an H/H2 transition 
layer below which the gas is molecular. The H2 v=l-0 S(l) 
line is emitted from the H/H 2 transition layer, which is close to 
the disk surface dStorzer & Hoilenbachlll999l) . See Fig. 13 in 
iDullemond et alj d2007l) for an illustration of the stru cture of 
a photoevaporative disk. IStorzer & Hollenbachl d 1999b pointed 
out that the ionization structure in a PDR with photoevaporation 
may differ from a classical PDR because the material is not at 
rest. They also found that the H 2 v=l-0 S(l) and v=2-l S(l) 
lines are mainly collisionally excited, and that fluorescence 
contributes only a small amount to the line intensities. The 
/v=i-os(i)/A=2-iS(i) ratio is typically 5-10. 

If the molecular hydrogen in the H/H2 transition layer gets 
hot enough to overcome the escape velocity it may contribute to 
the gas flow away from the disk. This could explain the small 
blueshift of the line profile in Fig. [4] 

5.2.2. Mass loss 

It is interesting to consider the scale of the mass loss due to pho- 
toevaporation. A charac teristic radius for thermal evaporation is 
dDullemond et al.ll2007l) : 



(5) 



where the sound speed equals the escape speed. However, sim- 
ulations have shown that evaporation is initiated at smaller radii 
than this analytical result sugge sts. That is, photoev aporation 
happens outside r cr > 0A5r g dDullemond et alj 120071) . Using 
M=2M and T=10 4 K as a typical temperature of the ionized gas, 
we find r cr ~ 3AU or about 20 mas. It thus seems very likely that 
the H2 emission caused by irradiation is accompanied by a pho- 
toevaporative wind although it is not directly detected in these 

ob servations. 

iHollenbach etalJdl994h derived an estimate for the hydro- 
gen mass loss rate due to photoevaporation, 



hid 



4 x 10" 



10 41 s- 



Me Msyl " 



(6) 



where d) is the stellar ionizing flux. Alexan der et al.l {2005) 
found O ~ 10 41 - 10 44 photons s _1 in five CTTs with ages 
of around 10 6 years. Assuming that the ionizing flux of T Tau 
N also falls within this range, we estimate the mass loss rate 
to be 5 x 10~ 10 - 2 x 10~ 8 Moyr _1 . The mass loss rate from 
a photoevaporative wind would thus be significantly smaller 
than the accretion rates inferred from SED modeling (Table [2]), 
while the upper limit is si milar to the obser ved accretion rate of 
(1 .4 - 5.9) x 10~ 8 M o yr _ 1 dBeck et al.ll2004l) . Since the accretion 
rate from the SED might be overestimated (Sect. [H the photoe- 
vaporative wind may have a significant effect on the dispersal of 
the T Tau N disk at its present evolutionary state. 

5.3. Other indicators 

Two scenarios seem likely to explain the H2 emission: shocks 
from a stellar wide-angle wind interacting with a flared disk, or 
irradiation from the central star onto the disk. The latter may 
or may not be accompanied by photoevaporative mass loss. The 



simplest way to distinguish between shock and PDR excitation 
is that the latter will create an atomic ionized layer giving rise to 
hydrogen recombination lines. Furthermore, a PDR will display 
strong emission from [OI] 63/im, 145/im and [CII] 158/jm as 
well as [Fell] 1.26/mi, 1.64/zm and [OI] 6300A and [SII] 6730A 
in high density regions. We examined the spatial distribution of 
the continuum-subtracted Bry line to search for extended Bry 
emission in comparison to the presumably pointlike PSF of the 
adjacent continuum. We did not find any conclusive evidence 
for any spatial exten sion of the Bry line. This is consistent with 
iKasper et al.l d2002i). who found that extended Bry emission, if 
any, is confined to within 6AU of the star. This argues against 
the existence of an ex tended PDR. 

On the other hand. lvan den Ancker et a li d 1999b found strong 
[OI], [CII] and [Fell] lines which they were unable to fit with 
shock models. Although these data suffer from low spatial res- 
olution (> 20"), these findings s uggest the pr e sence of a PDR 
component. With 2"resolution ISolf & Bohml d 19991) found a 
compact emission region in the optical [OI] and [SII] lines which 
is centered close to T Tau N and is characterized by near-zero 
radial velocities. Our SINFONI observations of the T Tau sys- 
tem included J-band data which are not presented here (see 
Gustafsson et al, in prep). These measurements are centered on 
T Tau S and do not cover T Tau N. Nevertheless, the data show 
indirect evidence of strong [Fell] 1.26/vmand [OI] 1.13/im emis- 
sion close to or originating from T Tau N itself. In the J-band 
data, continuum emission from T Tau N is scattered into the field 
of view and we find [OI] and [Fell] emission with the same spa- 
tial distribution as the scattered light. We therefore believe that 
the emission originates in T Tau N . The [OI] line traces the ion- 
ization front dMarconi et al.l!998h and the presence of these lines 
indicates that a PDR contribution may be present after all. 

One possible solution to reconcile these apparently contra- 
dictory observations is that a PDR exists but is confined to within 
~ 6AU from T Tau N. The infrared H2 emission analysed here 
extends to much larger radii and would then be due to shocks. 
In any case, the emission seems to be linked to a nearly face-on 
disk. Future observations with high spectral and spatial resolu- 
tion will be needed to confirm or disprove these suggestions. 
Such measurements should include deeper observations of the 
infrared H2 lines to constrain the line ratios, spectrally resolved 
measurements of the velocity, a search for extended hydrogen 
recombination lines, and observations to spatially constrain the 
[OI], [Fell], [SII] and [CII] lines. 

6. Implications for T Tau S 

T Tau S is located at a projected distance of 0Y7 (~ 100 AU) 
south of T Tau N. The extent of the H2 emission, as well as 
the SED modeling, indicate that the outer radius of the disk 
around T Tau N is 85-100 AU. Disks around stars in binary 
or multiple stellar systems will be truncated by their mutual 
gravitational influence. The size of the disk depends on the 
mass ratio of the stars and their separation. The disk size typ- 
ically ranges from 0.3-0.4 times the separation for mass ratios 
of 1-0.3, even when the disk and orbital plane are n ot coplanar 
( Artvmow icz & LubowHl99"4l: lLarwood et al.lll996l) . If the true 
size of the disk around T Tau N is ~ 100AU, the argumentation 
above implies that the orbital distance to T Tau S is ~ 300AU 
and that the inclination of the orbit to the line-of-sight is ~ 20°. 

If, on the one hand, T Tau S is located i n front of T Tau N, UV 
irradiation from the 2.7M star T Tau Sa dDuchene et alj |2006) 
could in principle contribute significantly to the H2 excitation 
on the side of the T Tau N disk facing toward us. However, in 
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that case the H2 emission would be accompanied by extended 
Bry emission as discussed in Sect. 15.31 This is not seen. On the 
other hand, if T Tau S is located behind T Tau N, the size of the 
disk indicates that it may very well contribute to the extinction 
of T Tau S. Given the small angular separation of T Tau Sa-Sb 
of ~0f.' 1 the disk around T Tau N would most likely obscure both 
Sa and Sb by an equal amount of material. Therefore, it does not 
help to explain why Sb is a normal T Tauri star and Sa show s 
the characteristics of an infrared companion dDvck et al.| [T982). 
The relative position of T Tau N and S and their outflows are the 
subjects of a subsequent paper (Gustafsson et al, in prep). 

7. Conclusions 

We detect emission from the H2 v=l-0 S(l) rovibrational line at 
2.12/im in a ring-like structure very close to T Tau N. We find 
that the weak H2 emission is most likely linked to a nearly face- 
on flared disk. Another possible solution is that the H2 emission 
originates from shocks impacting on the lower density walls of 
an envelope cavity. This scenario, however, requires that the high 
velocity jet in T Tau N be less collimated than in other T Tauri 
stars. The radius of the disk is ~ 85 - 100AU, based on SED 
modeling and the extent of the H2 emission. The velocity in the 
vicinity of T Tau N is consistent with the rest velocity of the 
star to within the errors. Both shocks associated with a wide- 
angle wind impinging on the disk and UV + X-ray irradiation 
from the central star onto the disk are plausible excitation mech- 
anisms which can reproduce the H2 flux. Both these mechanisms 
require a substantial disk around T Tau N. However, models and 
observations indicate that irradiation from the central star can- 
not excite H2 at radii much larger than 20 AU. Thus, the most 
likely excitation mechanism of H2 is that of a wide-angle wind 
impinging on a flared disk. A PDR created by irradiation may 
exist within ~ 6AU from T Tau N. 
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